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Abstract 

Two-leg t — J ladders are investigated in the framework of a combination of 
the phase string formulation and bond-operator representation. We develope 
a mean-field theory in the strong rung interaction regime, i.e. J± S> J, t, 
which provides a unified description of the undoped insulating phase and the 
low doping phase — the so-called CTS'O phase. Both of them are character- 
ized by the resonating-valence-bond (RVB) order parameter, with gap opened 
up in all spin excitations. The ground state of the doped phase is intrinsi- 
cally a superconductor with a d-wave symmetry, which is driven by the RVB 
correlations. The ground-state energy is in good agreement with numerical 
results. Phase separation is shown to occur beyond some critical value of J/t 
for given doping concentration. We also show that the spin gap in the doped 
phase is determined by quasi-particle-like excitations. The local structure of 
hole pairs as well as the spectra of various spin and charge modes are analyzed 
in comparison with other approaches. 
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I. INTRODUCTION 



Since the discovery of high-T c superconductors, the study of t— J type models has become 
an important topic of strongly correlated electron systems. The t — J model provides us one 
of the simplest examples of the nontrivial interplay between charge and antiferromagnetism 
in a doped Mott insulator. It is a problem generally difficult to tackle by conventional 
many-body methods due to the nature of strong correlations. 

Recently, ladder systems have received intensive studies both theoretically [|T|-|T3H and 
experimentally |H-ff7l. From theoretical point of view, the ladder t — J systems may be 
easier to investigate both numerically and analytically than the two-dimensional (2D) case 
related to the high-T c cuprates. However, the former may already catch some key physics of 
the latter and offer some important insights into the competition between charge and spin 
correlations beyond one- dimensional (ID) geometry. 

Like in ID and 2D, the physics of ladder t — J systems has been more or less well- 
understood at half-filling where only the spin degrees of freedom are present. For 
example, in the two-leg ladder problem, the ground state may be visualized as a spin liquid 
state of the condensate of short-ranged RVB spin singlets . In the strong rung interaction 
limit, a bond-operator representation provides a very useful description of such a 
spin liquid state. 

Also similar to ID and 2D, the doped case in ladder t — J systems poses a real challenge 
to analytic approaches: The central issue is how to correctly handle the competing charge 
and spin correlations once holes are introduced. In the two-leg ladder system, it has been 
established on numerical basis that on the small doping side (5 < 0.5) there exists a so-called 
C1S0 phase [§,|§ where spin excitations are all gapped (denoted by 50) while the density 
fluctuations of holes represent the only gapless mode there (denoted by CI). Such a phase 
is a superconducting phase with a d- wave- like symmetry. With the increase of the ratio J/ 1 
(where t and J are two parameters of the t — J model to be defined later), eventually a 
phase separation is found beyond some critical value of J/t. 
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Previously, a mean-field theory of lightly doped two-leg ladders proposed by Sigrist et 
al. |7]] showed that the short-ranged RVB state evolves into a superconductor with modified 
d-wave symmetry. It also gave a continuous evolution of the spin gap with doping. However, 
this theory fails to capture many features of the ground state and the excitation spectrum 
uncovered later by numerical work P~5|] . The main reason is that the fermionic RVB 
mean-field theory employed in this approach is not a very good description even at half- 
filling. It is the purpose of this paper to develop an analytic framework that can give a 
unified and systematic description of the aforementioned physical properties in both the 
undoped and doped phases of the two-leg ladder t — J model. We shall start from the strong 
rung interaction limit where it is natural to adopt the previously mentioned bond-operator 
representation [@,@ for the undoped case. The corresponding mean-field theory gives 



a reasonably good description of the ground state and spin excitations at half-filling as 
memtioned before. In the doped case, additional bond operators are apparently needed and 
they can be classified as the rung hole pair (d) and quasiparticles (a CT and a a ) besides the 
spin singlet (s) and spin triplet (t a ) operators. They compose of a complete basis convenient 
for describing the low doping case. 

For any dimensionality, it is highly nontrivial to get access to the doped phase from the 



undoped insulator of the t — J model. This is due to the fact that the Marshall signs fl9 
hidden in the half-filled spin background will be generally "disordered" by the motion of 
holes, leading to the phase string effect f2"0fl . Such a phase string effect cannot be repaired 
through spin flip processes as the latter always respects the Marshall sign at low energy. 
It implies that the nonrepairable "phase strings" left on the hole paths will be present in 
the ground state of the doped case. These "phase strings" play a role similar to the Fermi- 
surface phase-shifts originally proposed by Anderson. Indeed, in ID case, the phase string 
effect leads to the Luttinger liquid behavior. The phase string formulation |20[ provides a 
systematical method which reproduces correct exponents of various correlation functions. 



Nontrivial phase string effect in 2D mean-field theory has been also investigated in Ref. PT 



Incorporating the phase string effect thus becomes a necessary step to construct a sensible 



mean-field theory in the study of the two-leg ladder systems. 

Starting from the half-filling where the mean-field theory is based on a RVB characteri- 
zation with an order parameter (sj) = s, we are able to generalize the theory to the doped 
regime after incorporating the phase string effect. We find that the ground state in the 
doped phase is naturally a superconductor with a d-wave-like symmetry, as the consequence 
of the RVB correlations in the insulating phase. The mean-field ground-state energy is in 
good agreement with the numerical one at the doping concentration 5 < 0.5. Moreover, 
an instability of phase separation occurs in our mean-field state as the ratio J/t increases 
beyond some critical value, also consistent with numerical results. The present mean-field 
theory thus accommodates the most important physical properties of the doped two-leg 
t — J model previously identified only numerically. We would like to point out that without 
explicitly dealing with the nonlocal phase string effect at the starting point, a mean-field 
treatment would lead to a phase which is always unstable against phase separation, similar 



to the spiral instability in 2D case [22 



Furthermore, a series of detailed features obtained in various numerical work are also 
reproduced at this mean-field level. The energy spectra of magnons and quasiparticles have 
been determined in the C1S0 phase where they all exhibit finite gaps varying with the 
doping concentration. The gap of the former continuously evolves from the insulating phase 
while the one of the latter arises from the formation of Cooper pairs between quasiparticles. 
The minimum gap of creating a pair of quasiparticle excitations is smaller than that of the 
magnon. This indicates that the spin gap in the two-leg ladder system is generally associated 
with quasiparticles and shows a discontinuous evolution with doping, as first pointed out by 
Troyer et al. || . We also examine the local structure of hole pairs and show that the pairing 
on diagonal sites occurs simultaneously with the condensation of rung hole pairs. This point 



was also noted previously by Sierra et al. [11]. Since the pairing between quasiparticles 



also results from the condensation of rung hole pairs, a relationship between the pairing on 
diagonal sites and the spin gap in the two-leg ladders is then established. 

The rest of the paper is organized as follows: In Section 2, we introduce the phase 



string formulation and bond-operator representation. The general feature of the resulting 
Hamiltonian is discussed and the mean-field treatment is presented in Section 3. In Section 
4, we present our numerical analysis of the mean- field equations. The finial section is devoted 
to a conclusive discussion. For the sake of compactness, some useful and relevant formulae 
are listed in the Appendices. 

II. MATHEMATICAL FORMULATION 

We start with the original t — J Hamiltonian 

R t -j = P s {-tEfe + H.c.) + J]T(S, ■ S, - l -n inj )}P s (1) 

(*.i> (*.i> 

where P s is the projection operator which imposes the no-double-occupancy constraint such 
that the electron occupation number rii < 1. The conventional way to handle the constraint 
is to introduce the so-called slave-particle representation of electron operator: c JiCr — > f^bj^ 
such that the constraint rij < 1 is replaced by an equality condition: fjfc + Y^aHaha = 1- 
For our purpose, in the following we will use the slave-fermion representation in which fj and 
6j j(T satisfy the canonical anti-commutation and commutation relations, respectively. Then 
the t — J Hamiltonian H t -j = H t + Hj can be rewritten as follows: 

Hj = ~E b W-aV^- (2) 

z (id) 

In obtaining the above expressions, a replacement was made: bj j(7 — > (— aybj^ where 
a = 1,-1 for spin up and down, respectively. In this representation, the matrix element 
of Hj always remains negative-definite which is equivalent to say that the Marshall sign 
|T9| has been built into the basis |[20|| . But then the sign a appearing in H t indicates that 
holes dislike the Marshall sign hidden in the spin background, and their motion generally 
creates Marshall-sign mismatches on their paths known as phase strings. Since Hj respects 
the Marshall sign rule at low energy, the phase strings cannot be repaired through the spin 
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flip processes. Such a phase-string-type doping effect has been argued p(J to be the key to 



understanding the evolution of the ground state at finite doping. In the following, we first 



give a brief review of the phase string formulation developed in Ref. j2(J to deal with this 
nonlocal singular phase effect. 

A. Phase string formulation 

The basic idea underlying the phase string formulation is to "gauge away" the original 
singular source of the phase string effect shown in H t [Eq. (|2|)] such that the resulting form 



of the Hamiltonian becomes treatable in a perturbative scheme. According to Ref. p0|| , this 
procedure can be realized through a unitary transformation: 



U exp {--5>&(Q(1 - n? - £ <)} (3) 



where nh and n h j a are the number densities of holons and spinons with spin a at site j. 
Under the unitary transformation (H), the electron operators become 

cj,„ - hfbjA-vy*"" (4) 

where 

^^^expl^^WlE^t- 1 )}- 
bj,a = 6 3 >exp{--cr£^(0nJ, 

hj = f. exp {-z^^K}. (5) 

^ is the total number of holes. It is easy to verify that hj is a hard-core boson, i.e. they 
satisfy the following commutation relations: [hi, hj] = = i ^ j and {hi, hi} = 0, 

{hi,hf} = 1. Even though in the original definition, 6j(l) = Im ln(zj — z{), the choice of 
9j(l) is equivalent to a kind of gauge fixing and we will explicitly write down our choice later. 
With eqs. (ffl) and (El), the t — J Hamiltonian becomes 
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H* = -tY,{{j A «>%{j aA ^^ic + H.c.}, 

(id) 

Hj = ~ EC^JbLbJ^Ce^)^^ (6) 



where the gauge phases A{j and A^- are defined by 



M, = \nm)-oM (e< 



4 = *Eft(O-0*(OH fc - (7) 



In the case of one chain, we can choose 9i(l) such that all these gauge phases vanish |20 



Thus, all important phases are absorbed into eq. (f|) in a form of phase shifts [20] . This is the 
prefect case that the phase string effect is completely "gauged away" from the Hamiltonian. 
We shall see later that the situation slightly changes in the ladder case. 

We now focus ourselves on the two-leg ladder and define 9j tm (l,n) as follows: 

9j, m (hn) = 0, j > /, 

7T, j < I, 
7T 

Qj,m(j,n) = -, m = l,n = 2, 

-~ m = 2,n = l (8) 

where m, n G {1,2} are labels of legs. 1 and 2 indicate the upper and lower chains, 
respectively. This convention fixes the gauge phases as 

A ) h j+i = 4 ( n i,2 + n i+i,2)> 

A % + i = li 1 ~ S h ~ S hi,i) (9) 
where vh m and 5? are the number density of holes and the z- component spin operator 



at site j and chain m, respectively. represents the gauge phase on chain m. In 

the derivation of eq. we have used two identities: expi7r(l — 2S'? m + cn^ m ) = 1 and 
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expiir(a — <j')n h j m = 1, which are valid in the physical Hilbert space. With the help of eq. 



3,m 



f), we can rewrite the t — J Hamiltonian on two- leg ladders as follows: 



= E(£&A-a)EWW) 

3 ® <r' 
— £ (^j,m,o-^j,m,o-^j+l,m,-(T^i+l,ni,-(T) 

2 Z^l e U j,l,a U 3,l, -o- W j + l,l,-o-°J + l, 1,(7 

< AnA*}- (10) 

Note that in eq. ([H]), the coupling between different legs is explicitly distinguished from 
the one on the same leg. In the following, we will reformulate it in the bond-operator 
representation under the implicit assumption that J± ^> t, J. 



In contrast to the results of the single chain (see Ref. pOfl), there are some phase factors 
left in the Hamiltonian. If we increase the number of chains to infinity, they will turn out 
to become topological gauge fields. It is these gauge fields that strongly affect the dynamics 
of holons and spinons in two dimensions . On ladders, there are more than one path to 
connect two points while there is only one way to do it on the single chain. Consequently, 
in general, there is no closed path in the latter case [23]. That is why we are unable to see 



those "gauge interactions" arising from the phase string explicitly present in the ID t — J 
model. As for the ladder case, there is no way to completely gauge away these phase factors 
no matter how we choose the 9j^ m (l,n). Nevertheless, those phase factors in eq. ( ]T0| ) will 
become trivial in the bond-operator representation introduced below. 
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B. Bond-operator description 



In half-filled case, the two-leg ladder t — J model at J_l ^> J can be well described || 
at the mean-field level based on the bond-operator representation originally introduced by 



Sachdev and Bhatt [18]. Such a description can be easily generalized to the doped case as 
emphasized in the Introduction, At each rung, the physical Hilbet space is spanned by nine 
states which can be generated by applying the bond operators to the vacuum state | <po) as 
follows: 

dj\(j)o) = hf 1 hf 2 \0), 



a .j 



;„ | O > = (-cry +1 hlb^ a | 0), a+ CT | O > = (-^'&+i>+2 I °>. 

4I^) = ^(IT!)+IIT)), 

Here | a, a') = b^^b^.a* I 0) an d | 4>o) is annihilated by these bond operators, s represents 
the spin singlet and t a with a = ±1,0 represent spin triplet excitations. a a and a a particles 
carry the same quantum numbers as electrons, d particles are spinless and charge two. We 
choose s and t operators to satisfy canonical commutation relations as in Ref. || while d, 
a a , and a a operators are hard-core bosons because they contain h operators. The prefactor, 
(— iy , in eq. (|TT|) arises from (— cr)-? in eq. The no-double-occupancy condition is 

replaced by the following one: 

s t 8 i + J2 ttotj,* + Y,( a t,* a j,v + S>Oj>) + djdj = 1. (12) 

a=±l,0 a 

By using eq. (|11|), one can express the bilinear operators composed of h and b a by 
these bond operators. Their detailed forms are left in the Appendix A. We now obtain a 
Hamiltonian in which those gauge phases can be evaluated explicitly by substituting eq. 
lT|) into eq. (0). The resulting Hamiltonian can be divided into three parts: 

H 4 _j = H + Hx + H 2 
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where 



H = -N(^ + 1) ■ iy,cr(a p; a, n + H.c.) 

j 

j 

•( a /+v a J> + s i+i,<r s j>) + H - c -} 

~T XX( a j>J> + S>^^)4+i^+i + 0' ^ i + !)} 
j 

j 



^jitfSj + *J a *J,a + a £r°i,ff + a J«r 4,«r + dj^j ~ 1) (13) 



and 



-9 E{(E ii^yE ^^) + ( a -> a )} 
-rEft^+A+i' ( 14 ) 

Here /i is the chemical potential of holes. A, is the Lagrangian multiplier to impose the 



constraint (|T2|). N is the number of sites for the single chain. The attraction between d 
particles comes from the — l/AniUj term while those among quasiparticles arise from the 
exchange term. All the remaining terms are collected in H2, which consists of terms that 
involve either spin-flip processes or the creation or annihilation of triplet excitations. We 
shall see later that all spin excitations are gapped. Thus, we expect that the inclusion of 
H2 is supposed not to change the main features of our results and will not consider it in the 
following calculations. We leave the detailed form of H2 in the Appendix B and define our 
working Hamiltonian as H = Ho + Hi in the following sections. 
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III. MEAN-FIELD THEORY 



A. The low energy effective Hamiltonian 

Before taking the mean-field approximation, we would like to re-organize the general 
form of H t _j in the phase-string and bond-operator representation [Eqs. and (O)]. 

Then some features of it can be more easily revealed. 

We first express the a-operators in terms of the bonding and anti-bonding operators as 
follows: 

where a± j(7 denote the bonding and anti-bonding operators, respectively. Since a± i(T operators 
are still hard-core bosons, one may introduce the following Jordan- Wigner transformation 
to transform them into fermions without changing the Hamiltonian: 

o+ a = ef a Uj, at >j>0 <. r: r r (16) 

Here Uj = exp {in J2i<j,a( a + 1 o a +,i,<j + a ^ j a a -,i,o)}- and e^ a become fermions and satisfy 
the canonical anti-commutation relations. 

Secondly, because under the unitary transformation (Q), St becomes 

S !a = ( 1)^;,A,.. ^p {m EKi + nja) - i-fifo, 

Kj Z 

= ( l) i+1 bf^b jM exp {m E«i + < 2 ) + (17) 

Kj / 

instead of (— l) J 6^6j^, the t operators defined based on b^ a in flllD do not form a vector. 
Therefore, the coefficients before the terms tj+^o^o and tj+i >cr tt a in eq. (0) have different 
signs. The manifest rotational symmetry can be easily recovered by introducing the following 
unitary transformation: 

f t - t ex p{- CTff EKi + <2)}- (is) 

Kj 

for tj >a with a — ±1. 
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Then H t _j is simplified after substituting eqs. ([15j) and ( p!6|) into eqs. (|i~3|) and (|T4] ) and 



performing the transformation (18) [25]: 



and 



Ho = + i - 1 E(4.^> - <;,/,.,) 



4 2' ; 



2 

j 

+(x + ^ 5>£^> + fe) + + j) E 44 

i j 

_ T 12{( e l* e j>° + ^j,cr)d'} +1 d j+1 + (j <-> j + 1)} 

j 

J 

Z 3 



E A i(4 s i + *Ja*j,a + 4<^> + + d+c^ - 1) (19) 



3, a, a' 

+(e+ +1>ff et_ CT - et + ve+_ ff )(e i ,_^e i+ i )ff i - e 3 - - CT >e i+ i, CT /)} 

(2°) 

z i 

In addition to the spin rotational symmetry, translation symmetry, and the electromag- 
netic U(l) symmetry, we note that the Hamiltonian of the two-leg ladder is also invariant 
under the exchange of chain indices to which we can assign a parity operator. Under the ex- 
change of chain indices, the bond operators transform as follows: s — > s, t a — > —t a , d — > d. 
For quasiparticles, it is a a <-> — aa a , or e a — ► e a , e a — ► — e CT . We see that magnons and 
quasiparticles in the anti-bonding band are parity odd while d bosons and quasiparticles in 



the bonding band are parity even [£6[ . 
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We shall see later that both quasiparticles and magnons have gaps at low doping con- 
centration. Consequently, we can integrate them out (treating s as a c- number.) and obtain 
the low energy theory described by the following hard-core boson (HCB) model: 

H = -** E(4+i rf i + H - c -) + V E d i d j d j+i d j+i ( 21 ) 
i i 

where t* is the effective hopping amplitude and the dominant contribution to V comes from 

the attraction between d particles in eq.(|20|). This model can be solved by a bosonization 

approach [^]. For V < —2 | t* | the system is phase separated and this occurs only at 

very large values of J± for physically reasonable values of J/t. For example, J± > 31. 8t 

for J/t = 0.5 (see Troyer et al. ||). In the region where the system is stable against the 

phase separation, the low-lying excitation will be the phase fluctuations of <i-particles, which 

corresponds to the collective charge mode. 

Although the HCB model, eq. (|21|), is appropriate to describe the low energy properties 

of two-leg ladders in the lightly doped region, it can not address questions such as the 

internal structure of hole pairs, which is related to the nature of the superconducting order 

parameter, and how the spectra of those gapped modes vary with the hole concentration. 

Later we will answer these questions by a mean-field treatment of eqs. ( |T§D and 

B. Mean- field equations 

To proceed with the mean-field approximation, we first note that the undoped two- 
leg ladder is characterized as a spin liquid with non- vanishing RVB order parameter (sj). 
Following Gopalan et al. ||, we take the ansatz for the spin part as: (sj) = s, (tj, a ) = 
0, Qa = (tj+i,att a ) 5 an d Pa = (ttatj,a) where a = ±1,0. Because of the rotational 
symmetry, Q + = Q_ = Q and P + = P = P . We set P = P + + P_ + P and Q = 

Q+ + Q- + Qo- 

For the charge part at finite doping, we define the following mean-field parameters: 
Xa = ( e j h +i, - e i,o-) an d Xo- = (ej+i,crZj,o-) • Again due to the rotational symmetry, x+ — X- 
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and x+ — X-- We define x — X+ + X- an d X — X+ + X-- Moreover, we take Xj = A in 
accordance with the translational invariance along the chain direction. 

Note that there is a linear <i-operator appearing in the first term of Hi [Eq. (pOpj. 
describing the process that a rung hole pair dissolves into two quasiparticles or a pair of 
quasiparticles is recombined into a rung hole pair. Such a term (with Sj — > s) looks similar 
to a Hubbard- Stratonovich transformed four-fermion interaction in the BCS theory with d 
playing the role of the order parameter. Thus, we will treat dj as a c-number by assuming 
dj = d and neglect its phase fluctuations. A solution with non-vanishing d will then im- 
mediately lead to forming Cooper pairs for e-p and e\ according to the first term in Hi. To 
make things more transparent, let us define the following pairing fields for quasiparticles: 

Then 

A + A = - E( a i+i,-4-a - 

A ~ A = ~ 12 a ( a i+l^t-a ~ ^t+l,a a t-a)- 

Here A + A and A — A represent the hole pairing along the chain direction and diagonal 
sites, respectively and are shown schematically in Fig. 1. Now it is easy to see that the 
pairing on the rung and diagonal sites must occur simultaneously. This is what the first 
term of Hi tells us and it is consistent with the picture emerged from numerical studies 
0,0 ■ 

Naively, it seems that A and A are not tensors under the spin rotation and thus our 
mean-field ansatz may break the rotational symmetry. This is, in fact, disguised by the phase 
string effect as we have discussed before. Our mean-field ansatz indeed respects rotational 
symmetry which is left to be discussed in the Appendix C. 

Based on the above mean-field ansatz, we finally obtain the following mean-field Hamil- 
tonian: 

n MF = N{\-±- J -- (p ± + \)s 2 + ( X + x)Qt + ^(A + A) 2 

14 



~ix - X? - + 2a* + A - ^ - J)d 2 } 
-(*(<? + ! ^) - ^(x - X)) E(4+i,^> + ^- c -) 

3 

j 

+(-t - /i - A + -f + 3 - - J -d 2 ) £ e+ a e ha 

i 

3 

-d(A + A) + V2tsd) + H - c -) 

3 

-(^(A + A) - V2tsd) YM+xA.- + H - c -) 

3 

+ 1 E((^ 2 " t(x + x)KJi+i,* + Js 2 tj,Jt + i,- a + H.c.} 

Z 3 

+(^"A)E^- (23) 

3 

Eq. ( p3j ) can be diagonalized by Bogolioubov transformations. We leave the procedure in 
the Appendix D and write down the diagonalized Hamiltonian in the following: 

R MF = fl + EW^ oc k + PtPk) + E k {ata k + f%fa)} 

k 

+E (24) 

k 

where 

= E(^fe - E k - E k ) + - 2// - ^ + 3 - + ( X + x)Qt ~ {\j± + A)s 2 

+^(A + A) 2 - Jq(x - X? ~ + 2/i + A - (25) 

is the mean-field ground state energy. (In fact, Qq is the zero temperature grand potential.) 

The parameters d, s, //, A, x> Xi A, A, ^> an d Q are determined by solving eq. (p2|), the 
following self-consistent equations: 



Q — E^'+^ia)' P ~ E(*i«*i«)' 

x = E( e jWv)' x = E(4+i^)' ( 26 ) 
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and the saddle-point equations: 



dtt 



0. 



0. 



dQ 
ds 

dn 



0. 

-2N8, 



(27) 



dd ' dfx 

where 5 is the hole concentration. We obtain the following mean-field equations from eqs. 



, (ID, and (|3): 



1 - s 2 + d 2 
l + dP-5 



P + 25, 

ilV(-^ + ^), 
2 N , E k E k 



(28) 
(29) 



(2 i (x + x)+2 ^_£_ 1 + ,)J 
^ + )-Q + A(, + ,) + ^(A-A) 



-^-S(A-A), 
N 



(30) 
(31) 



A 

X 

Q 



l 

iV 



i ^ . . r* 



AT 



sin 



Ei. 



- —7 / COSK — , Y = — — > COSK-=- 



31^ , A fc 



P 



2N^ u k 



3 
2' 



(32) 
(33) 
(34) 



Notice that eq. ( P0| ) says that A — A = as long as d — 0. 

We close this section by a remark. When 5 = 0, i.e., the undoped case, our mean- field 
equations are not exactly reduced to those in Ref. |J. The difference arises from the fact 
that there are three components for t operators and each contributes J2k to VLq. Thus, 
the zero-point energy is Y,k f^fc instead of J2k \^k- A similar effect also appears in eq. (|25|) 
in which the coefficient of NX changes from 3/2 to 1/2. This point was neglected in Ref. 
||. We will see later that this difference increases the spin gap in comparison with that 
obtained in Ref. [Hi. 



IV. RESULTS 
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A. Undoped case 



As a reference point, we first examine the results of our mean-field equations for the 
undoped case. By defining A k = A k | x=0 , the relevant integrals in the mean-field equations 



can be expressed as elliptic integrals |28 

1 



N 



E 



A fc 



1 



1 

N 



cos k 



k U k 

Ak - 2n fc 



{- 



7T VI + V 



2v 



v 



1 2v_ 



)}• 



KU 



2u 



' 2u 
1 + u 



)} 



(35) 



where K(£) and E(£) are respectively the complete elliptic integrals of the first and second 
kind with modulus £. The dimensionless parameter v is defined in the Appendix D. 

Then, in the undoped limit, eqs. (p8j), (p9|), (p0[), fl3ll) , (|32|) , (0), and (0) are reduced 
to the following forms: 



* - * = i-{ -J_ «r (, 



+v / l + ^( 



' 2i/ 
1 + v 



)}, 



(36) 



3 _A_ 

4 Jj_ 7TZ/ Vl + ^" 'V 1 + ^' 

-vT+i>£( 



2i/ 



)} 



(37) 



where r] — J/J±. The spin-triplet excitation spectrum is given by 



1 A 



OJk — J ±(-7 — )Vl + ^COsfc. 

4 J i 



(38) 



In eq. (|38|), we have to assume < z/ < 1; otherwise the mean-field equations would break 
down. This is verified in the following calculations. The band minimum is at k = it and the 
spin gap is determined by 



V. 



(39) 



17 



We can analytically study the asymptotic behavior of the spin gap for small values of r] 
(and hence v) to see the difference between our equations and those in Ref. |J. For small 
values of £ the elliptic integrals K(^) and E(£) can be expanded in a power series as 

K(0 = f(l + i? 2 + ^ 4 --), 

and we obtain 

v = 2,(1 + -irf + Ofe 4 )), 

and the spin gap is given by 

A t = J ± (l-r ] + ^ V 2 + 0(r ] 3 )). (40) 

We find that our result, eq. (POD, is larger than that in Ref. and is much closer to the 
strung rung interaction result, which is Jj_(l — r\ + |r/ 2 + 0(r/ 3 )) [|J. At the end of Sec. Ill 
we have pointed that there is a numerical factor missing in the zero-point energy in Ref. f| 
which is responsible for this discrepancy. 

To obtain A t at any value of J/ J±, we numerically solved eqs. (p6|) and fl37|) . The results 
are shown in Fig. 2. The spin gap at the instropic point, r] — 1, is about 0.501J, which is 
very close to the numerical result — 0.504 J 0. Of course, the mean-field approximation 
is only justified in the strong rung interaction regime and we do not expect the theory to 
be extended into the region with r] > 1 where the coupling between spins on the same leg 
becomes dominant over the rung coupling. In fact, our calculation shows that the spin gap 
continuously increases beyond rj > 1, but according to Ref. [3], it should smoothly diminish 
to zero as r/ approaches 0. 

B. Phases in doped case 
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1. The CI SO phase 



The phase diagram obtained from numerical studies P|3]T^] shows that at most values 
of J/t the two- leg ladders fall into the universality classes of Luther- Emery and Luttinger 
liquids for small and large doping concentration, respectively. The former and the latter are 
respectively denoted by C1S0 and C1S1 phases. (CmSn means that there are m gapless 
charge modes and n gapless spin modes pUfl .) 

The present theory using the bond-operator description presumably works at small dop- 
ing for the Luther-Emery type phase. A mean-field solution with non-vanishing d is found at 
5 < 0.5 which is stable against the phase separation when J/t is not too large. We interpret 
this mean-field state as the C1S0 phase. This is because in this situation there are gaps in 
the spectra of both quasiparticles and magnons. This can be understood as the following: 



To let quasiparticles be gapless, T k and T k must be zero f29|| . From eq. ([32|), it results in 
vanishing A and A. This implies that d = due to eq. (|30D . Therefore, a non-zero d will 
induce a gap for quasiparticles. The gap in the spectrum of magnons is a continuation of 
the one in the undoped case, which ensures the validity of the extension of the mean-field 
ansatz from the undoped case to the finite doping. The only gapless excitation in this region 
is the density fluctuations of hole pairs as discussed before. 

Ground state energy: To examine the validity of our mean-field ansatz, we first 
compute the ground state energy and compare it with numerical results. The mean-field 
ground state energy per site, E , is given as the following: 

3 

= + »s. (4i) 

In the above derivation, we have used eq. Also note that the number of sites is 2N. 

We calculated E with J/t = 0.5 and its doping dependence at various J±/J's is plotted in 
Fig. 3. By comparing with the numerical results in Ref. |TI ], we find that both the tendency 



of the energy versus the doping concentration and its absolute magnitude agree well with 
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the data obtained by the recurrent variational ansatz (RVA) as well as the density matrix 



renormalization group (DMRG) method. For comparison, in Fig. 3 the DMRG results [II 



at 5 = 1/8 and 1/2 are also shown. (Note that we set J = 1 in Fig. 3 while J = 0.5 in Ref. 



IIJ but J/t = 0.5 is the same.) We see that the agreement is especially good in the region 
with large values of J±/J and small doping concentration as expected for the bond-operator 
representation. (The comparisons with the RVA results are even better over the whole 
5 < 0.5 region.) Such a good agreement over a wide range of parameters indicates that 
our mean-field treatment based on the phase string and bond-operator formalism indeed 
captures the basic physics of doped two-leg ladders in the strong rung interaction regime. 
In the following, we focus on some detailed properties by solving the mean-field equations 
at J/t = 0.5 and J ± /J = 10. 

Local structure of hole pairs: Next, we would like to discuss the local structure of 



hole pairs in the CI SO phase. As has been discussed in Ref. |llj, holes will form pairs along 
the diagonal sites as well as along the rung and chain directions. This diagonal pairing is 
energetically favored by the t term and the most probable configuration of two dynamical 
holes in a two- leg ladder [[II]]. In our formalism, the amplitudes for pairing along the rung, 
diagonal, and chain directions can be respectively represented by d, A — A and A + A. We 
calculated (| A — A \)/d and (| A + A \)/d and plot them in Fig. 4. We found that both 
decrease with increasing 5. In the low doping region, the amplitude of hole pairs on diagonal 
sites is almost comparable to the one of rung hole pairs and always larger than that on the 
chain direction. The amplitude of hole pairs along the leg being smaller than other hole 
configurations reflects the fact that the rung bonds are stronger than the leg bonds in the 



underlying two-leg spin ladder. The above results were also pointed out by Sierra et al. fllT 
They proposed a dimer hard-core boson (DHCB) model to describe the low energy properties 
of two-leg ladders in the strong rung interaction regime, which contains both the charge and 
spin degrees of freedom in contrast to the HCB model. The bond-operator formulation also 
retains these high energy modes. In particular, the coexistence of the diagonal pairing and 
rung hole pairs is further manifested in our approach. 
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Pairing symmetry: We also calculated the expectation values of pairing fields along 
rung and chain directions. The corresponding operators are defined as follows: 

1 



A x (j) = ■j=Y^ (TC i+W>« c 3&-<» 



\ti) = -j= I] vc jX<J c j)2 ,-o. (42) 

These pairing fields can be represented by bond operators and we list their explicit forms in 
Appendix E. In terms of the above mean-field parameters, their vacuum expectation values 
are given as follows: 

i{A x ) = L(A + A)[s 2 -2d 2 -^-- 7 

t td 
+Q- 7 (X + X)-V / 2^(A-A)], 
J Js 

i(A v ) = ±ds. (43) 

The results are shown in Fig. 5. It is clear that there is a critical hole concentration 5 C . 
(In our case, 5 C = 0.37 for J/t = 0.5 and J±/J = 10.) In the low doping regime 5 < 5 C , 
the pairing symmetry shows d-wave-like behavior while for 5 > 5 C , it becomes s-wave-like 



symmetry. The difference can be attributed to different internal structures of hole pairs |TT 
When 5 < 5 C , holes doped into a spin liquid state with RVB correlations form pairs with 
d x 2_ y 2-Y\ke structure. However, in the overdoped region, one moves into the low density limit 
characterized by electrons doped into a background with an internal s-wave-like symmetry. 

Spin excitations: There are two kinds of excitations which carry non-trival spin quan- 
tum numbers. One is the magnon, which is represented by t operators in our formulation 
and is the spin triplet excitation around q = (it, it). The other type of spin excitations are 
quasiparticles, which carry spin-1/2. The band minima of quasiparticles in bonding and 
anti-bonding bands are at q = (0,0) and q = (0,7r), respectively. The behaviors of their 
gaps varying with hole concentration are shown in Fig. 6. We found that the gap of magnons 
increases while quasiparticle gaps decrease with increasing 5. In addition, it is easy to see 
that the low-lying spin modes with odd and even parity are magnons and quasiparticles in 
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the bonding band, respectively. Real spin-1 excitations in two-leg ladders are composed of 
magnons or pairs of quasiparticles. The gap of the latter is still smaller than that of the 
former. Thus, the spin gap in the C1S0 phase is determined by quasiparticles instead of 
magnons. This was also shown by numerical studies ||. 



2. Phase separation 

The above-discussed (7150 superconducting phase may become unstable against phase 
separation when the value of J/t becomes large in the t — J ladders 0. The reason is that 
in the large J limit the gain in exchange energy by maximizing the number of AF bonds 
outweighs the cost in kinetic energy. 

The stability of the C1S0 solution against the phase separation can be examined by 
studying the compressibility k: 

k- 1 - 5 2 ^ 
~ d 85- 

At Jj_/ J = 10, we found that k diverges at J/t = 1.2 and becomes negative when J/t > 1.2. 
This implies that our mean-field solutions with uniform hole density become unstable against 
phase separation when J/t > 1.2 [see Fig. 7(a)]. For J±/J = 5, the situation is similar as 
shown in Fig. 7(b). The only difference is that the boundary between the C1S0 phase 
and the phase separation region is moved to J/t = 1.6. At J±/J = 2, the critical value 
of J/t where the phase separation occurs not only increases but also becomes strongly 
doping-dependent as shown in Fig. 7(c). The latter trend is quite similar to that found in 
numerical studies P,[IT| for the isotropic case, i.e. J±/J = 1. We note, however, that further 



reducing J±/J towards the isotropic limit in our mean- field theory does not improve more 
of the comparison with the numerical results since the ground-state energy starts to visibly 
deviate from numerical data at J±/J<2 even for small 5 as shown in Fig. 3. 
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V. DISCUSSIONS 



In this paper, we proposed a mean-field description on doped two-leg ladders in the 
strong rung interaction limit based on the phase string formulation. The transformation to 
bond operators is a natural choice in this formalism. With the help of bond operators, we 
can easily separate the high erergy and low energy processes in the t — J Hamiltonian. Thus, 
a mean-field treatment becomes straightforward. Naively, there are two competing pairing 
channels in eq. - the pairing between e CT and e_ CT and the one between e a and e_ CT (or 
e a and e_ CT ). However, the condensation of the hard-core boson d and a non- vanishing RVB 
order parameter s demand that the latter dominates. Furthermore, this type of pairing 
implies the formation of spin singlet bonds not only along the chain direction but also on 
diagonal sites. The latter turns out to be an important low energy structure in various types 
of t — J models ||10|| . This singlet bond becomes a strong nearest-neighbor one after one of 
the holes hops next to the other. As a consequence, the formation of this kind of singlets 
can maximize the hopping overalp with other hole configurations and lower its energy. This 
feature is a necessary result in our formulae as shown in eq. fl30D or the linear d term in 
Hi [Eq. (pODI- We have to emphasize that this mechanism for pairing comes from the t 
term and is quite different from the "broken-bond" effect though the latter does enhance 
hole pairing somewhat. Also, the pairing must cause a gap opened up in the quasiparticle 
spectrum. Thus, all spin excitations are gapped in this region. 

Here we would like to make some comments on the effects of phase string in two-leg 
ladders. If we directly apply the bond-operator representation to the original t — J Hamil- 
tonian without explicitly taking into account the nonlocal effect of phase string, then we 
would obtain a Hamiltonian with a different form in the i-term. For example, the t-term 
with a linear d in Hi would involve quasiparticle pairing between different bands which is 
a high energy process now. Subsequently, the pairing between quasiparticles would only 
come from the four-fermion attraction in the J-term of the Hamiltonian. If we still use the 
similar mean-field ansatz to treat this Hamiltonian, the results would be incorrect because 
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the compressibility is always negative. This implies that the solution is thermodynamically 
unstable. Such an instability is actually similar to the spiral instability in 2D case when 
one tries to generalize the Schwinger-boson mean-field theory to the doped case without con- 



sidering the phase string effect |?T] • Therefore, it is necessary to take into account the phase 
string effect in order to acquire a correct mean-field theory on the doped antiferromagnets 
regardless of dimensionality. 

Muller and Rice have suggested a possible C2S2 phase existing between the Nogaoka 
phase with very small vaules of J/t and the C1S0 phase with intermediate values of J/t 
\\12\\ . They provided some numerical evidence to support this conjecture. Physically, the 
appearence of this phase can be understood as the following: When the holons move fast, 
i.e. t 3> J, J±, the gain in kinetic energy may outweights the cost by breaking the rung 
hole pairs and rung singlets. Thus, the phase coherence between bonding and anti-bonding 
bands is lost and the two-leg ladder is effectively decoupled into two chains at low energy. 
If this picture is correct, then the bond operators (especially the d-operator) are no longer 
a good description of the low energy degrees of freedom in this region. On the other hand, 
as we have seen in eq. (P^), there are always attractions between quasiparticles on the 
nearest-neighbor sites arising from breaking the singlet bonds. If they are not completely 
compensated by some repulsive forces at least at the intermediate scale, the underlying 
magnetic structure may still be a gapped spin liquid with a small spin gap and the observed 
C2S2 phase perhaps is a finite size effect. But the present mean-field theory, which works 
in the limit J± 3> J, t, cannot be directly applied to this regime to address those issues. 

For large doping concentration, the anti-bonding band of electrons is empty and the 
system falls into the C1S1 phase. To describe this phase, the bond-operator representation 
is not convenient. We have to go back to eq. ([H]). Nevertheless, the problem that there are 
non-trivial phase factors in the Hamiltonian rears its head again. These phase factors are 
the interactions arising from the phase string effect and entail careful treatment. Otherwise, 
important physics may be lost. The pursuit along this direction is beyond the scope of the 
present paper. 
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APPENDIX A: THE BOND-OPERATOR REPRESENTATIONS OF BILINEAR 

OPERATORS IN THE HAMILTONIAN 

In this Appendix, we list the bond-operator representations of those bilinear operators 
appearing in the Hamiltonian in the following: 

(-a) j+1 

h] )X b 3Xa = -/=— a+_ a ( Sj + <Tt jt0 ) + (-ay +1 a+ a t j!a + (-ayd+a j>a , 

(-a) j 

hf !2 bj,2,a = —j=-al_ a (sj - (Ttj, Q ) + (-a) J a+ a t j!a + {-a) 3+1 d^ a j>a , 
2&j,i,A2 -<r = (-l) 3 y/2sj, 

a 

Z a=±l,0 a'=±l 



~' t ~ a j,a a j,o-i 

2 ! 



a=±l,0 a'=±l 



-\-(lj a Q,j^ c , 

(-1V 
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+ (-l) i+1 aJ ff a i> _ ff . (Al) 
The above equations should be considered only in the physical Hilbert space. 

APPENDIX B: THE FORM OF H 2 

We list H2 in the following: 

H 2 = 77 E a ( s j+i t t,o + t i+i,o s j h )( a j h +i, (T a 7> - 

E^+l.O^ + dj+ltjfi)(aj' + l,a^t~^ + tij+l^Oj-o) 

-itJ2(tj+i,adj - d j+1 t j:a )(a+ +l a a+ a - a+ +l a a+ a ) 
3 

~i E a {( s i*Jo + *j,o4)( a i+i,«r a i+i.» ~ + (i ^ i + !)} 

i 

+ ^7 E CT {( E V)( a j+i,- a i+^ + S + +i A+i,<r) + U J + !)} 

s j CT '=±i 

X! (^,0^+1,0^^+1 -a ~~ ^,o4+l,o4«^'+ 1 ,a) 
Z j,a=±l,0 

-(J ^J + l)} 

+tE( E E tf+i^W) + #- c - ( B1 ) 

4 j cr=±l cr'=±l 

APPENDIX C: A AND A ARE SPIN SINGLETS 

We would like to discuss the transformation properties of A and A under the spin rotation 
and verify that they are spin singlets. Due to the same reason as that for t operators, a a 
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and da denned in ([TT|) are not spinors. Therefore, we are unable to directly conclude that A 
and A are not tensors under the spin rotation. To make the rotational symmetry manifest, 
we perform the following unitary transformation: 

4 Z Kj 

at - 4. ex P {^7 " *l EKi + < 2 )>- (CI) 
Under the above transformation (|CT|), A + A and A — A become 



A + A -f i J] <r(aj+ir a t-° ~ a t+i,<r a l-*) > 

cr 

a - a - e + ,)• (C2) 

cr 

Now it is clear that A and A are indeed spin singlets. We also examine other terms in H 
with eq. ( JC1| ) and confirm that our mean-field ansatz respects the rotational symmetry. 

APPENDIX D: DIAGONALIZATION OF EQ.(|23|) 



Here we present the procedure to diagonalize eq. (|23|) . After performing Fourier trans- 
formations on all operators by Oj = ^= J2k Oke tkXj , the mean- field Hamiltonian in eq.(p3[) 
becomes 

R MF = N{\- J -±- J -- (| J ± + \)s 2 + ( X + x)Qt + ^-(A + A) 2 
-{- Q (x-x) 2 -(^ + 2» + \-±-J) ( P} 

k 

+ ^2{iT k e^et ki + iY k e^et kl + H.c.) 
k 

+ T,i A kti a tka + n k (tiA_ a + **«*-*-«)} (di) 



where 



e k = -{t{2a* + s 2 + Q)- - x)) cask 
-t-fi-X + ^ + ^l-d 2 ), 
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t k = -{t{2d 2 + s 2 + Q) + -(x - X)) cosk 

T k = (j(A + A) + 2V2tsd) sink, 

f fc = (—(A + A) — 2V2tsd) sink, 

A k = (Js 2 - t( X + X)) cos k + ^- - \, 

= — s cosfc. 
2 

Here the lattice spacing has been taken to be unity. 

Eq. ( pi|) can be diagonalized by the following Bogolioubov transformations: 

otk = u k e k ^ - v k et kl , /3 k = u k e- kl + v k e^ p (D2) 
afc = u k e k1 - v k e± kl , (3 k = u k e- k[ + v k e kp (D3) 
7fc Q = cosh 9 k t ka + sinh d k t\_ a . (D4) 

The coefficients u k , v k , u k , v k , cosh^, and sinh^ are given by 

u k = cos<p k e 1 ^, v k = sin0 fc e"^, 

u k = cos0 fe e l 4, v k = sin0 fe e _l 4, 

cos 2 fc = + ^-), sin 2 (p k = i(f - ^-), 
i hj k z hj k 

cos 2 <p k = -(1 + sin 2 fc = -(1 - ^-), 
z hj k A £j k 

E k = y/? k + VI E k = sjel + Tl (D5) 
cosh 2 ^ = i(^L + 1), sinh 2 ^, = - f), 

Z LO k Z LO k 

u k = J A 2 - 4IT 2 . (D6) 



In eq. ( |JJ6| ) , both A k and A 2 , — 4IT 2 , have to be positive. This constraint will be enforced in 
our numerical analysis. If we define v = (x + X)t/(^i — ^) an d v = 2 Js 2 / (^4 — A), then the 
band minimum of t particles occurs at k = n when v < u/2 and at k = when v > vj2. 
With the help of eqs. (|D2"1) , (p3|) , and ( |D4D , we obtain eq. ([24]) . 
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APPENDIX E: THE BOND-OPERATOR REPRESENTATIONS OF PAIRING 

FIELDS 

Here we give the bond-operator representations of pairing fields in the phase string 
formulae. First, we need the corresponding representation of electron operators. They are 
as the following: 

7T 7T 

Cj,i, a = exp H-(l + a)) exp {i- ^ (2S£ m - 1 - an l>m )} 

l>j,m 

{^=a^__ a {sj + at jiQ ) - ia^tj^ + d+a jta }a Nh , 
c jt 2,a = exp 0'-(l + a)) exp {i- J2 ( 2S i,m ~ 1 - <™j,m)} 

l>j,m 

{-^pf+_ a {s 3 - t7t A0 ) + />/•,,/,„ + d+a^}^. (El) 
By pluging the above formulae into eq. (0), we obtain 

l>j,m 

A.(j)=exp(2« 2 S^J 

l>j+l,m 

1 + 

~~^J2^3 tj+l,v®'j+l,a a j-o- + ^•+l^i,-o-%+l,o-0 J) _ (7 ) 

+-o-[dJ fl a i+ i i0 .at (T (s i + o"f ii0 ) - d+a+ +1(T a ji(T (s i+ i + o-^+i.o)] 
+-aj+i iCT «j^ J >('Sj+i + o"*j+i,o) - tj+\,<r{sj + &tj,o)] 
+7^tf +lia d+_ a (s J+1 + o-t i+ i )0 )(s j - (T^-o)}. (E2) 



We have to emphasize that it is necessary to take into account a Nh in eq. (|E1|). Otherwise, 
the representations of pairing fields would be incorrect. 
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symmetry forbids the condensation of d particles. However, d can condense after taking 
into account the phase string. This will result in very different physics. 
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Figure Captions: 

Fig. 1 The local structure of hole pairs: (a) diagonal site, (b) chain direction. The solid 
line and open circles represent the spin singlet bonds and holes, respectively. 

Fig. 2 The spin gap, A t , as a function of 77 = J / J± where E = A t / J±. 

Fig. 3 The ground state energy per site with J/t = 0.5. Different curves correspond to 
J±/ J = 1, 2, 4, 6, 8, and 10. The data denoted by the cross are obtained with DMRG 



and taken from Ref. 11 . 



Fig. 4 The weight of different types of hole pairs with J/t = 0.5 and J±/J = 10. The solid 
line and open circles represent the weights of hole pairs on diagonal sites and chain 
direction, respectively. 

Fig. 5 The vacuum expectation values of pairing fields as functions of the doping concen- 
tration at J/t = 0.5 and J±/J = 10. The dashed line and open circles represent A x 
and A y , respectively. 

Fig. 6 The gaps of spin excitations at J/t = 0.5 and J±/J = 10. The solid line, open 
circles, and stars correspond to quasiparticles in bonding band, anti-bonding band, 
and magnons, respectively. 

Fig. 7 The boundary between the CIS'O phase and the phase separation region with (a) 
J ± /J = 10, (b) J ± /J = 5, and (c) J ± /J = 2. 
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